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ABSTRACT 

Among various analytic approximations for the growth of density fluctuations in the expand- 
ing Universe, Zel'dovich approximation and its extensions in Lagrangian scheme are known to 
be accurate even in mildly non-linear regime. The aim of this paper is to investigate the reason 
why these Zel'dovich-type approximations work accurately beyond the linear regime from the 
following two points of view: (1) Dimensionality of the system and (2) the Lagrangian scheme 
on which the Zel'dovich approximation is grounded. In order to examine the dimensionality, 
we introduce a model with spheroidal mass distribution. In order to examine the Lagrangian 
scheme, we introduce the Fade approximation in Eulerian scheme. We clarify which of these 
aspects supports the unusual accuracy of the Zel'dovich-type approximations. We also give an 
implication for more accurate approximation method beyond the Zel'dovich-type approxima- 
tions. 

Subject headings: cosmology: theory — galaxies: clustering — gravitation — large-scale 
structure of universe 



1 Introduction 



As an analytic approximation for the growth of density fluctuations in the expanding Universe, 
Zel'dovich approximation (ZA hereafter) is known to be unusually accurate even in mildly 
non-linear regime for unknown reason (Zel'dovich 1970; 1973). Recently the extensions of this 
ZA (Post- and Post-post- Zel'dovich Approximations, PZA and PPZA hereafter) have been 
developed (e.g., Bernardeau 1994; Buchcrt 1994; Bouchet et al. 1995). These Zel'dovich- 
type approximations are confirmed to be far better than any other analytic approximations 
in Munshi, Sahni k Starobinsky (1994), Sahni & Shandarin (1995) and Sahni & Coles (1996) 
using a spherical-collapse model. However so far, there has been no clear explanation why 
these Zel'dovich- type approximations are so accurate beyond their range of applicability. In 
this paper, we would like to address this issue from limited aspects and with limited models. 

In the Zel'dovich-type approximations, we work on the Lagrangian coordinate scheme in 
which the location of a mass element x of the fluid is expressed by the initial location q and 
the time dependent displacement vector ^ as 

xiq,t) = q + ^{q,t). (1.1) 

Then the density contrast S[x{q, t),t] = det[dxi/dqj]~^—l is determined by solving the equation 
of motion 

Vj* = 2ff2!!{(a;,t), (1.3) 

where V^: is the spatial derivative with respect to Eulerian coordinates x. This nonlinear 
equation for ^ can be solved by the method of iteration considering d^i/dgj to be a small 
parameter. The first iteration corresponds to the ZA. With further iterations, PZA and PPZA 
are considered to improve the accuracy. 

There seem to be at least two possible grounds for the Zel'dovich-type approximations: 

(1) In the plain parallel mass distribution, equations (1.2) and (1.3), with S = 1/(1 + — 1, 
reduce to 

(-) 

where = d^/dq. This differential equation (1.4) is exactly the same as that for the 
linear density contrast Si,, thus ^' oc (5l. Therefore, in this case, ZA becomes exact at 



2 



least before shell crossings occur. This one-dimensional-exact property is considered to 
support the validity of this approximation also in three-dimensional systems. 

(2) ZA is unique in the sense that it is based on the Lagrange coordinate scheme while all 
the other nonlinear approximations are based on the Eulerian coordinate scheme. This 
fact may make the ZA extraordinary excellent. 

In this paper, we would like to clarify the reason why Zel'dovich-type approximations are so 
good from the above limited points of view. 

In order to elucidate the aspect (1), we introduce a model of spheroidal collapse. In the 
previous analytical work (Munshi et al. 1994; Sahni & Shandarin 1995; Bouchet et al. 1995) 
only spherical symmetric models, which are analytically solvable, have been used to examine 
the validity of the Zel'dovich-type approximations. By changing the axes ratio in our spheroidal 
model, we can freely control the effective dimensionality of the system. For example, a pro- 
late collapse is effectively dimension two and an oblate collapse is effectively dimension one. 
Therefore, if the one-dimensional-exact property in the first aspect (1) in the above gives the 
very reason for the Zel'dovich-type approximations, their accuracy would be the best in oblate 
collapse and may be better in prolate collapse than in spherical collapse. 

In order to explore the aspect (2), we introduce the Fade approximation in Eulerian co- 
ordinate scheme. This is not a simple polynomial expansion with a small parameter but a 
rational polynomial expansion. In this paper, we compare this sophisticated approximation in 
the Eulerian scheme with Zel'dovich-type approximations in Lagrangian scheme. If the aspect 
(2) gives the very reasoning for the Zel'dovich-type approximations, their accuracy is far better 
than this Fade approximation. 

In the course of our study, we need to establish the hierarchy of accuracy in various nonlinear 
approximations. Therefore we examine the other known nonlinear approximations in Eulerian 
scheme such as (a) linear perturbation and higher order perturbation approximations, (b) 
frozen-flow approximation, and (c) linear-potential approximation, as well as Zel'dovich-type 
approximations and the Fade approximation. 

The organization of this paper is as follows: In section two, we summarize various nonlinear 
approximations in gravitational instability theory. They are applied to the spheroidal pertur- 
bation model in section three. The first aspect is examined in this section. In section four, we 
introduce the Fade approximation in Eulerian scheme. The second aspect is examined in this 
section. In section five, we clarify which aspects are supported from our analysis. In the last 
section six, we conclude our analysis and briefly mention the possibility of the approximation 
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scheme beyond the Zel'dovich-type approximations. 



2 Nonlinear Approximations in Gravitational Instability Theory 
2.1 Equations of motion and linear perturbation scheme 

In the gravitational instabihty theory, the non-relativistic matter with zero pressure in an 

Einstein-de Sitter universe is described by the following set of equations (see e.g., Peebles 
1980), 

(5 + V ■ [(1 + (5)v] = 0, (2.1) 

V + 2Hv + {v ■ \/)v + V$ = 0, (2.2) 

= ^H^5, (2.3) 

where x, v{x,t), ^{x,t) are respectively position, peculiar velocity, peculiar potential in comov- 
ing coordinate, which correspond to ax, av, in physical coordinate. An over dot denotes 
the time derivative and V = d/dx denotes the spatial derivative with respect to comoving co- 
ordinates. The scale factor a varies as a oc t^/^ and the Hubble parameter is H = a/a = 2/{3t). 
Although we consider Einstein-de Sitter universe throughout this paper, it is straightforward 
to generalize our analysis to J7 7^ 1, A 7^ universes. 

For these non-linear equations, we show various approximations in the following. In the 
linear regime {S ^ 1), we can safely neglect the nonlinear terms and obtain relatively simple 
solution (see, e.g., Peebles 1980). Neglecting decaying mode, the solution is 

dj^{x,t) = ^Si^{x), (2.4) 

$L(cc,i) = ^^'A-IJl, (2.5) 
VL{x,t) = -^V^L, (2.6) 

where is the inverse Laplacian: 

These quantities simply evolve as Si, oc a, $l o"^, vl oc a~^/^. 

For the later convenience, we introduce the peculiar potential <p = a^$. In Einstein-de Sitter 
universe, the linear peculiar potential is constant, 0L(a;,t) = const. = (l){x). 
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2.2 Higher order perturbation methods in Eulerian scheme 



In higher order perturbation methods, nonlinear correction terms are added to the linear so- 
lution; S = 6l + (5^^) + 5^^) + • • •, ^> = $L + + ^^^^ + • • •, V = vi^ + + + • • •, 
where (^f") <I>("^ is assumed to be of order This expansion enables us to solve the 

nonlinear equations (2.1)-(2.3) order by order. The second-order expression of the perturbation 
is relatively simple (e.g., see Peebles 1980; Fry 1984): 

5 2 
S = Si, + -6l + S^iifm + -(pL,ij(Ph,ij (2.8) 

where (p = A~^Si,. However, the expression of higher order is complicated. For the detailed 
expression of the third and fourth-order solution in Fourier space, see Goroff et al. (1986). The 
explicit expression for the higher order is not necessary for our analysis, so is not quoted here. 

2.3 Frozen flow and linear potential approximations 

Matarrese et al. (1992) introduced the frozen flow approximation (FF, hereafter). In this 
approximation, the velocity field v{x, t) is kept fixed to the value of linear perturbation scheme: 

VFF {x,t) = -^^"^Mx), (2.9) 

A particle in FF moves simply along the line determined by the fixed linear velocity fields (2.9), 
i.e., the position of each particle, x{t), is described by the differential equation, 

where d/dt is the Lagrangian time derivative. It is convenient to rewrite the equations as 
follows: 

fa = -S^T^^'*^*^'- 

On the other hand, Brainerd, Scherrer &; Villumsen (1993) and Bagla & Padmanabham 
(1994) introduced linear potential approximation (LP, hereafter) which is based on the assump- 
tion that the gravitational potential evolves according to the linear perturbation scheme. As 
a result, particles effectively move along the lines of force of the initial potential ^in. Thus the 
position of each particle x{t) is described by the differential equation, 

d'^x ^^dx 1 „ , , , 

This equation is also rewritten as follows: 

d'^x 3 dx 1 „ , , , „ , „x 
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2.4 Zel'dovich approximation and higher order perturbation methods in Lagrangian 
scheme 



In the Lagrangian perturbation methods (see, e.g., Bernardeau 1994; Buchert 1994; Bouchet 
et al. 1995; Catclan 1995) we consider the motion of mass elements labelled by unperturbed 
Lagrangian coordinates q. The comoving Eulerian position of mass element q at time t is 
denoted by x{q,t). The displacement field ^{q,t) defined by equation (1.1) is the dynamical 
variable in this formulation. Taking divergence and rotation of the equations of motion (1.2), 
(1.3), we obtain the equations of motion for ^(q,t), 



dt 



dt^ 



2H 



d^. 



3,1 



dt 



V J Ik 







(2.14) 
(2.15) 



where d/dt is Lagrangian time derivative, Jij = dxi/dqj = Sij + W^ij, J = det Jij, and we used 
the relation (Vx)i = {J~^)ji{^q)j- In usual treatment of Lagrangian perturbation methods, one 
assumes an additional condition, i.e., vorticity-free condition Vx x v = 0, which is equivalent 
to 

which replaces equation (2.15). The solutions with this vorticity-free condition form a subclass 
of the all general solutions [rotational perturbation is argued by Buchert (1992) and Buchert & 
Ehlers (1993)]. Therefore, equations (2.14) and (2.16) are solved perturbatively for derivatives 
of displacement field = ^^^j '^^if ' ' ' keeping only terms of leading time-dependence 

[sec Buchert & Ehlers (1993), Buchert (1994), Bouchet et al. (1995), and Munshi et al. (1994) 
for detail]. In Einstein-de Sitter universe, the time dependence of each terms is separated from 
its spatial dependence: 

/ 2 \" 



The resulting perturbative solutions, up to third-order, are 



(1) _ 



(2) 



-di(t>Liq), 

14 



(3) 



(2.17) 

(2.18) 
(2.19) 

(2.20) 
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The first-order solution is equivalent to the ZA. While the first- and second-order displacement 
field is irrotational, the third-order solution has both the rotational part (the last term) and 
the irrotational part (first two terms). 

3 Spheroidal perturbation 

In this section, we examine our aspect (1), onc-dimensional exact property of Zcl'dovich-type 
approximations applying various non-linear approximations summarized in the previous section 
to the spheroidal collapse model. 



3.1 Equation of motion of an ellipsoid 

A particle motion in an Einstein-de Sitter universe is described by the following equation of 
motion: 

+ 2H— = -V^^x,t), (3.1) 



dt"^ dt 

V^^ = ^H'6{x,t), (3.2) 

where d/dt is the Lagrangian time derivative as before. In a homogeneous ellipsoid, the density 
perturbation 6{x,t) is given by 

^(,,<).4(*)e(i-^-^-^), (3.3) 

where ctj are the half-length of the principal axes of the ellipsoid and 9 is a step function. 
The solution of the Poisson equation (3.2) inside the homogeneous ellipsoid is known (sec, e.g., 
Kellogg 1953; Binney & Tremaine 1987), 

A-iG fl-4-4-4)=7E ^i^l (3-4) 
{r.h.s. is only for inside the ellipsoid) 

where 



roc "J 

Ai = aia2a3 / (af + A)"^ J| (a/ + X)~^''^dX. 



(3.5) 



These coefficients Ai automatically satisfy the following constraint 



3 

Y,Ai = 2. (3.6) 

i=l 
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Thus, the solution of equation (3.2) becomes 



^ = ^HX^A{t)x^. (3.7) 

i=l 

The quadratic form of the potential (3.7) implies that there exists a solution that retain 
homogeneity of the ellipsoid if the initial velocity field is linear in space and if the homogeneity 
of outer region is assumed. (Lynden-Bell 1962; 1964; Lin, Mestel &; Shu 1965; Icke 1973) 

Combining the equations (3.1) and (3.7),we obtain 

Xi + 2Hxi = —H^d^AiXi, (3.8) 

which describes the motion of the particles inside the ellipsoid. Generally, the variables Se 
and Ai^s depend on the position of other particles, i.e., volume and shape of the ellipsoid. 
Fortunately in our model, it is sufficient to consider only three particles at the coordinates 
{ai{t),0,0), (0, a2(t),0) and (0,0, a3(t)) because these three particles completely characterize 
the motion of the entire ellipsoid. The density contrast of the ellipsoid is given by 

(5e = 4(1 + 5i,)^Hil^Hi^ - 1, (3.9) 

where ain = Q;(iin), Oin = o,{tm) and tin is the initial time. The equation of motion for the three 
points are given by 

ai + 2Hai = - ^H^eAiai. (3.10) 

The above equations (3.5), (3.9) and (3.10) are the closed set of equations of motion which 
describe the motion of the entire ellipsoid. Because of the simplicity, the homogeneous ellipsoid 
model can be used as an approximation for protoobjects in structure formation in the universe 
(White & Silk 1979; Eisenstein & Loeb 1995; Bond & Myers 1996). 

We numerically solve these equations. For the initial condition of the numerical integration, 
we adopt Sin = 10~^-10~^ and the ZA for velocity field. For sufficiently small value of Sin, it 
is accurate enough to adopt ZA to prepare the initial velocity field. Actually, adopting PZA 
or PPZA does not change the results and the simple ZA is sufficient. In this paper, we adopt 
spheroidal symmetry, ai = a2, in which case, AiS have analytic forms as follows (Peebles 
1980): 

A, = A2 = ^{l + h), As = ^{l-2h), (3.11) 
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where 



h 



■ sm e — 



2e2 ' 



e=\ 1 



as 



3-e2 



e= \ 1- 



(ai =a2> as) 



(ai =a2 < as) 



(3.12) 



This spheroidal symmetric model is sufficient to study the dimensionality of the system. 



3.2 Linear perturbation scheme 

The evolution of spheroidal perturbations in linear perturbation scheme is simple: 



which reduces to 



5j^{x,t) = —5i^ e 1 



Si,{x,t) = ±a eil 



a- 



inl 



a. 



in2 



a 



in3 , 



t2 



X2 
2_ 



t2 



(3.13) 



(3.14) 



Hill "in2 "inS / 

after the normalization of scale factor Qj as fljji — l^inl* Here and after, upper sign corresponds to 
the positive perturbations 5c > 0, and the lower sign corresponds to the negative perturbations 
5c < 0. Solving Poisson equation by equation (3.4), linear peculiar potential inside the spheroid 
is given by 



(1 + hin)ixi + Xi) + (1 - 2hi^)xi 



(3.15) 



In Figures 1-5, we plot Si, against the numerically solved true density contrast Stme for various 
axis-ratios and for positive and negative perturbations. 



3.3 Frozen flow approximation 

We now derive the evolution of spheroidal perturbation in FF. Linear peculiar potential inside 
the spheroid is given by equation (3.15) and the differential equation (2.11) reduces to 

^^-^-¥^(1 + M^(i,2), § = ^^(1-2^X3. (3.16) 



da ' 3' ' da 
where (1,2) denotes the subscript 1 or 2 in this order. The solution of these equations are 

1 - 2h 



X{i,2) = 9(1,2) exp I =F — - — a 



X3 = q3 exp =F 



(3.17) 
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where g^'s are constants of integration which correspond to the initial condition, xi qi as 
a — 0. Because 5 — as a — 0, and p oc (a;iX2X3)~^, evolution of the density contrast is 



1 = exp(zba) — 1. 



(3.18) 



It should be noted that the solution (3.18) does not depend on the parameter hi^, and has the 
same evolution as in the spherical model, = 0. In Figures 1-5, we plot the evolution of FF 
with thin dot-dotted dash lines against (^true- 



3.4 Linear potential approximation 

Similarly, the solution of LP for homogeneous spheroids is obtained as follows. Equations (2.13) 
and (3.15) imply the equation of motion of particles in LP: 



d X(i 2\ 3 dx(i 1 , 



da2 

d^xs ^ 3 dx^ 



da? 2a da 2a 
The solutions of these equations are 



(1 - 2/ii„)a;3 = 0. 



X3 

Xil,2) 
X3 



(1,2) 



v/2(l + /ii„)a 

93 



sin ^2(1 + hin)a 



= sin \/2(l — 2hi„)a 

^2(1 - 2/iin)a V ^ 

, sinh ^2(1 + h,^)a 
V2(l + /iin)a V 

93 



^2(1 - 2/iin)a 
The evolution of density contrast is, therefore, 

2^2(1 + /?,;„)(! -2/lin)l/2^'^/2 



sinh J 2(1 — 2h\n)a 



{S > 0), 



{S < 0). 



s= < 



sil? a/2(1 + /iin)asin a/2(1 - 2/iin)a 
2^/2(1 + /iin)(l-2/iin)V^a3/2 
. sinh^ -y/2(l + /iin)asinh -^2(1 - 2/iin)a 

In the spherical model {hin = 0), the expression reduces to 

\ 3 



1 (<5>0) 
1 i6<0) 



' (^) -1 (*>0) 
ysm v2a J 



2a 



-1 iS<0) 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



\sinhv2o^ 

which corresponds to the result already derived by Brainerd et al. (1993). In Figures 1-5, we 
plot the evolution of LP with thin short-dashed lines against 5true- 
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3.5 Zel'dovich approximation and higher order Lagrangian perturbation methods 



In ZA, the particle motion is described by the equation (2.18). Using equation (3.15) for 
spheroidal case, we obtain 



4'- 



T- 



9(1,2) , 



1 + hi, 



T 7i (1 - 2/iii 



(3.24) 



As for the second-order term (2.19), the inverse Laplacian can be calculated by equation 
(3.4), because we assume that there is no fluctuation outside the spheroid all the time, i.e., 
6 = outside the spheroid in q-space. Thus in the Lagrangian perturbation scheme, we can 
easily identify the boundary of the spheroid. This is the point that the Lagrangian perturbation 
scheme is technically advantageous than the Eulerian perturbation scheme. We will return to 
this argument later. Thus equation (2.19) implies 

V-JS) = -^^^^^^(1 + ^in - C - ^in), 



(3.25) 



28 



'-{l-2hin-ht^ + 2h] 



Similarly, we can calculate the third-order terms according to the equation (2.20): 



9(1,2) 



504 



(23 + 23hu 



S9ht„ - 25ht + 44/1? + SOh? 



=F 



504 



(23 - 46/iin - 39/if„ + 92/if„ + 2hi - 60/if^) 



(3.26) 



The density contrast in these Zel'dovich-type approximations is given by 

1 

1, 



where W'i = iZ^f ^ for ZA, = <Z^f ^ + ^ for PZA, and = <ff ^ + l^f ^ + ^^-f ^ for PPZA. The 
relation between ip and ^ is given by (2.17) 

In the spherical case, h-m = 0, the equation (3.27) reduces to 

-3 



^-1 



-1, 



a _ 23a^ \ 

^ 3 ~ 21 ^ 1701 J 



-3 



(3.28) 
(3.29) 

(3.30) 
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for, respectively, ZA, PZA and PPZA, which correspond to the result by Munshi et al. (1994). 
The Lagrangian perturbation methods for ellipsoidal collapse with respect to mass function is 
considered by Monaco (1997). 

In Figures 1-5, we plot the evolution of ZA, PZA and PPZA against Jtrue- 



3.6 Second- and third-order perturbation methods in Eulerian scheme 

In contrast to Lagrangian perturbation methods, the surface of the spheroid cannot be ex- 
plicitly expressed in Eulerian perturbation methods. This problem makes the calculation of 
the spheroidal perturbation difficult in Eulerian perturbation methods. We circumvent this 
difficulty by transforming the expression already obtained in Lagrangian perturbation scheme 
to that in Eulerian perturbation scheme. In doing so, we compare the small expansion param- 
eter in both schemes. It is in Lagrangian perturbation scheme and the density contrast 
6 in Eulerian perturbation scheme. We notice these two parameters are the same order, i.e., 
6 ^ 0{^ij). We also notice the n-th order perturbative solution in Eulerian scheme, (5^"^ is of 
order {6lT and ^-JJ ~ ©(^-Jj )" in La grangian scheme. Therefore ^-"^ ~ (5^"\ 

Thus, to obtain the perturbative series up to n-th order of density contrast in Eulerian 
scheme, it is sufficient to expand the density contrast of n-th order in Lagrangian scheme 
by parameters ^^^j and re-express it in terms of 5l- In our case, ^'^"^ oc o", so Eulerian 
perturbative series can be simply obtained by expanding equation (3.27) in terms of expansion 
factor a. The result is 

17 . 4 , , [Ml U , _£ 3 _^,4^„3 



S = ±a+ { — + —hf^] ± { + hi hi /ifJ « • (3-31) 

V21 21 'V V567 189 '° 81 189 ^ ^ 

Although we adopt the above method in this paper, this is not the unique choice. Actually, 
we can straightforwardly obtain the Eulerian perturbation expansion independently from the 
Lagrangian perturbation scheme. However in the Eulerian scheme, our assumption " there 
is no fluctuation outside of the spheroid" becomes ambiguous. This is because, in Eulerian 
coordinate scheme, the location of the surface of the spheroid cannot be explicitly described 
as mentioned above. In Lagrangian scheme, we assumed that there is no fluctuation outside 



the evolved surface of the spheroid in section 3.5. In Eulerian space, on the other hand, it 



seems technically favorable to assume that there is no fluctuation outside the initial surface 
of the spheroid. This is because the inverse Laplacian (3.4) is easily solved for the latter 
assumption but is difficult for the former assumption. The sole ambiguity in our model comes 
from how to fix the fluctuation outside the spheroid. This fixing is necessary for us to obtain 
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the numerically exact solutions for the evolution of density contrast. For comparison, we quote 
the result adopting latter assumption in Eulerian space, no fluctuation outside the initial suffice 
of the spheroid: 

/17 4 .\ , /341 38 . 80 16,4 \ o 

S = ±a+l — + —h? ]a^±( + hi /i? + ht a^. (3.32) 

V2I 21 'V V567 105 '° 405 945 'V ^ ^ 

These two expressions above are very similar in the following two aspects. Firstly the 
spherical terms (i.e., terms independent of hin) are the same in both expressions. This is because 
for the spherical perturbation, the evolution of density contrast at a point is determined only by 
the inside of the spherical shell on which the point is located; the artificial fixing of fluctuations 
outside the shell does not affect the evolution of density contrast at the point. Secondly the 
two expressions agree with each other up to the second-order. This can be understood as 
follows. In linear perturbation scheme, the artificial fixing of fluctuation outside the spheroid 
is not necessary simply because S oc a and the perturbation is uniform in the spheroid. In 
second-order methods, the fluctuation outside the spheroid does not generally vanish however 
we force it to vanish artificially. This artificial fixing does not affect the evaluation of the 
density contrast inside of spheroid within the second-order perturbation scheme because the 
second-order density contrast is determined as the non-local functional of density contrast of 
linear perturbation scheme. The artificial fixing, however, affects the evaluation of the third 
and higher order calculations. Thus the artificial fixing of fluctuations results in the difference 
of third and higher order perturbation terms. 

In Figures 1-5, we plot the evolution of second- and third-order result (3.31 thin graph) 
against t^true- 

4 Pade approximation 

In this section, we examine the aspect (2), Lagrange scheme, for the validity of Zel'dovich- 
type approximations. The Zel'dovich-type approximations are unique in the sense that they 
are grounded on the Lagrangian coordinate scheme. It is advantageous to use this Lagrangian 
scheme because the inertia term is linearized in velocity. Is Lagrangian scheme the indispensable 
reason why Zel'dovich-type approximations work accurately beyond the linear regime? For the 
purpose of addressing this issue, we introduce now the Pade approximation (see, e.g.. Press et al. 
1992) in Eulerian coordinate scheme. This is an approximation for some unknown underlying 
function f{x) with rational function whose power series expansion agrees with a given power 
series of /(x) to the highest possible order. Pade approximation can even simulate the poles of 
the underlying function f{x) and is generally better than simple polynomial approximations. 
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Figure 1: Density contrast in various approximations for spherical density perturbations, 
tical axis is the density contrasts (^approx of various nonhnear approximations and the horizi 
axis is the true solution ^true in spherical collapse model. The true solution corresponds t( 
diagonal straight line. The distance to this line represents the accuracy of each approxima 
Fig. 1(a) is for of positive perturbation and Fig. 1(b) is for of negative perti 

tion. Some of lines in this figure is previously appeared in Munshi et al. (1994) and Sah 
Shandarin (1996). 
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true true 



Figure 2: Density contrast in various approximations for oblate density perturbations, 
same figure as in Fig. 1, but for oblate spheroidal collapse. Fig. 2(a) is the case with the i] 
axis-ratio ains/cKini = 0.8, and Fig. 2(b) is with ams/oiini = 0.3. 
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true true 



Figure 3: Density contrast in oblate spheroidal collapse of negative perturbation. Same f 
as Fig. 2 but for negative perturbations. Fig. 3(a) is the case with the initial axis- 
Qiins/cKini = 0.8, and Fig. 3(b) is with ains/cKini = 0.3. 
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Figure 4: Density contrast in prolate spheroidal collapse of positive perturbation. The i 
figure as in Fig. 2, but for prolate spheroidal collapse. Fig. 4(a) is the case with the h 
axis-ratio ains/oiini = 1-2, and Fig. 4(b) is with ams/cKini — 3. 
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true true 



Figure 5: Density contrast in prolate spheroidal collapse of negative perturbation. Same f 
as Fig. 3 but for prolate spheroidal collapse. Fig. 5(a) is the case with the initial axis- 
Qiins/cKini = 1-2, and Fig. 5(b) is with ains/cKini = 3. 
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Pade approximation for a given unknown function /(x) is expressed as the ratio oi 
polynomials. 

M 

R{x) ^ ^ , 

1 + ^ bkX^ 

k=l 

where and bk are constant coefficients. Suppose we already know the first M + N coeffic 
Cfc of a series expansion of the function /(x) around x = 0. 

M+N 
k=0 

Then the above coefficients ak and bk are determined so that the first M -\- N coefficients 
series expansion of R{x) agree with the coefficients Cfc.! This condition 

Co = ao, 

N 

bmCN-m+k = -CN+k, k = 1, - ■ ■ ,N 

m=l 

k 

^bynCk-m = ak, k = l,---,N 

m=0 

fixes the coefficients ak and bk- 

Now we use Pade approximation with the perturbative expansions in Eulerian coord 
scheme. The density contrast in the spheroidal model up to the third-order perturbati 
given by equation (3.31). The corresponding Pade approximation is given by 



^ 17 + 4/^2 / 214 110,2 104,4 \ 2 

1 T -a + ht H ht H ht a 

21 V3969 1323 81 1323 'V 

The second-order Pade approximation is obtained discarding the term in the denomii 
of the above equation. 

These approximations are already superposed on top of the previous approximatioi 
Figures 1-5. At a glance, the Pade procedure dramatically improves the accuracy of the i 
perturbation approximations. This accuracy is almost the same order as Zerdovich-typ( 
proximations. 

^Choices A'^ = MorA'' = M + l are usually adopted. We used A'' = M + 1 in this paper. 
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The limit of one-dimensional collapse is achieved by setting e ^ 1 in equation (3.12 
h[n —1- In this limit, the equation (4.3) reduces to the exact solution, S = ±a/(l : 
Therefore the Pade approximation in the present model also has the one-dimensional-( 
property as ZA. i 

5 Why are Zel'dovich-type approximations so good? 

In this section, we now discuss the main theme of this paper "why are Zel'dovich-type apj 
imations so good?" based on the preceding two sections. 

As we can observe in Figures 1-5, the curves of approximations are almost smooth an( 
accuracy of each approximation can be observed at any density contrast. Therefore, we d 
to measure the accuracy of each approximation with each initial axis-ratio by the qua 
^approx/^true Sit Stme = 4, Well withiu the nonlinear regime. This measure of accurac 
logarithmic scale, is plotted in Fig. 6, where the horizontal axis is the initial axis-rat 
spheroidal perturbations. In the similar way, we plot the measure of accuracy ^approx/^tr 
^true = —0.6 for negative density perturbations (voids) in Fig. 7. 

These figures 6 and 7 summarize all of our analysis. We observe from them the folio 
facts: 

(i) Zel'dovich-type approximations are far better than any Eulerian approximations (e^ 
the Pade approximation) in the spheroidal model as in the spherical model. Higher ( 
Zel'dovich-type approximations are definitely better than lower order Zel'dovich- 
approximations all the time. 

(ii) The hierarchy in accuracy of each approximation remains the same for all initial axis- 
of spheroidal positive as well as negative perturbations . 

(iii) The accuracy of Zel'dovich-type approximations becomes better in both prolate 
oblate initial conditions compared with the spherical symmetric perturbations. Mor( 
the accuracy in oblate initial conditions is better than that in prolate initial condit 
All the other Eulerian approximations except Pade approximation, on the other 1 
have exactly the opposite tendency; the accuracy of them becomes worse in both pr 
and oblate initial conditions. Moreover the accuracy in oblate initial conditions is \ 
than that in prolate initial conditions. 

■'"The Pade approximation based on (3.32) also has the one-dimensional-exact property and has almoi 
same behavior as (4.3). 
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(iv) Pade approximations are better than any other Eulerian approximations. Higher ( 
Pade approximations are definitely better than lower order Pade approximations al 
time. Pade approximations also have the one-dimensional-exact property. 

(v) The accuracy of Zel'dovich-type approximations and that of Pade approximation: 
almost the same. 

(vi) Zel'dovich-type approximations and the Pade approximations approach the exact soli 
from opposite directions. For example in Fig. 6, all the Pade approximations ovei 
mate the exact solution while all the Zel'dovich-type approximations underestimati 
exact solution. Similarly in Fig. 7, PPZA overestimate the exact solution whil( 
corresponding third-order Pade approximation underestimate the exact solution. 

Results (i) and (ii) make us confirm the excellence of the Zel'dovich-type approxima 
also in the spheroidal perturbations. We also confirm the consistency of the higher ( 
Zel'dovich-type approximations; higher order iteration always yields better accuracy. R 
(iii) supports our first aspect that the validity of the Zel'dovich-type approximations is grou 
on the one-dimensional-exact property of them. Actually, the Zel'dovich-type approxima 
gradually become much accurate when the system gradually deviates from the spherical 
metry. They are most accurate in oblate collapse, which is effectively dimension one. ' 
are second most accurate in prolate collapse, which is effectively dimension two. Result 
reminds us of the potentiality and the consistency of the Eulerian scheme. Result (v) 
proves our second aspect that the validity of Zel'dovich-type approximations is grounde 
their Lagrangian scheme. Result (vi) signifies that Zel'dovich-type approximations and 
approximations are definitely different scheme despite their similar behavior in Figs. 6 a 
and their similar appearance in their expansions [(3.30) and (4.3)]. 

Therefore our results support the aspect (2) that the validity of the Zerdovich-typ( 
proximations is due to the one-dimensional-exact property of them and disfavor the aspec 
that the validity is grounded on the Lagrangian scheme of them. 

6 Conclusions and Discussions 

In this paper, we have discussed the validity of the Zel'dovich-type approximations fron 
two aspects; (1) the dimensionality of the model and (2) the Lagrangian scheme on v 
the Zel'dovich-type approximations grounded. We introduced a model of spheroidal 
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Figure 6: Density contrast of various approximations in spheroidal collapse of positive pe 
bations. This figure summarizes all the positive perturbations with various initial axis-rj 
The horizontal axis represents the initial axis-ratio a-ms/ami and the vertical axis repre 
the ratio of density contrast 5approx/<^true evaluated at Stme = 4. Both the axis ratio and de 
ratio are in logarithmic scale in this figure. The right square panel is the magnification o 
part surrounded by the dotted line. 
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Figure 7: Density contrast of various approximations in spheroidal collapse of negative 
turbations. This is the same as Fig. 6 but for negative perturbations. The horizontal 
represents the initial axis-ratio ains/ctini and the vertical axis represents the ratio of de 
contrast (5approx/^true evaluated at Jtrue — —0.6. Both the axis ratio and density ratio a 
logarithmic scale in this figure. The 1st ~ 3rd Eulerian approximations are so bad that 
are outside of this figure. 
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perturbations and compared several nonlinear approximations as well as Zel'dovich-typc 

proximations. 

We found the following facts from the spheroidal models. The accuracy of Zel'dovich- 
approximations and Pade approximations becomes better in both prolate (effective dimei 
two) and oblate (effective dimension one) initial conditions compared with the spherical 
metric perturbations. Moreover the accuracy in oblate initial conditions is better than 
in prolate initial conditions. These results may show that the Zel'dovich-type approxima 
are more accurate in lower dimensional collapse. All the other Eulerian approximation 
the other hand, have exactly the opposite tendency; the accuracy of them becomes wor 
both prolate and oblate initial conditions. Moreover the accuracy in oblate initial condi 
is worse than that in prolate initial conditions. 

The above facts are consistent with the aspect (1) that the validity of the Zel'dovich- 
approximations is due to the one-dimensional-exact property of them. On the other hand 
above facts conflict with the aspect (2) that the validity is grounded on the Lagrangian scl 
of them. For the final confirmations of the aspect (1), we need further analysis on much \ 
class of models. 

The last fact (vi) in the previous section suggests a possibility to construct a better ap] 
imation beyond Zel'dovich-type approximations by applying Pade method. Since Zel'do 
type and Pade approximations are definitely different things, we may obtain new inform; 

when the Pade procedure is applied on the Zel'dovich-type approximations. The considera 
on this possibility will be affirmatively reported in our separate publications. 

For the actual analysis of observational quantities in the Universe, it is advantageoi 
use Eulerian scheme rather than Lagrangian scheme because the technique of classical 
theory is most naturally applied in the former scheme. In this context, Pade approxima 
in Eulerian scheme, as they are accurate as Zel'dovich-type approximations, may open 
perspective in the analysis of growth of density fluctuations in the expanding Universe 
course we should keep in mind that the validity of the Pade approximation has not yet 
fully established and therefore we should avoid blind applications of Pade method in cosmc 

We thank Masaaki Morita for helpful discussions on the dimensionality at the occasi( 

Tokyo Seminar at Mitaka National Observatory May 1997. We also thank Oki Nagahar 
many suggestions in numerical calculations. 
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